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Abstract 

We observe mixing between two-electron singlet and triplet states in a double quantum dot, 
caused by interactions with nuclear spins in the host semiconductor. This mixing is suppressed 
by applying a small magnetic field, or by increasing the interdot tunnel coupling and thereby the 
singlet-triplet splitting. Electron transport involving transitions between triplets and singlets in 
turn polarizes the nuclei, resulting in striking bistabilities. We extract from the fluctuating nuclear 
field a limitation on the time-averaged spin coherence time T£ of 25 ns. Control of the electron- 
nuclear interaction will therefore be crucial for the coherent manipulation of individual electron 
spins. 



A single electron confined in a GaAs quantum dot is often referred to as artificial hydro- 
gen. One important difference between natural and artificial hydrogen, however, is that in 
the first, the hyperfine interaction couples the electron to a single nucleus, while in artificial 
hydrogen the electron is coupled to about a million Ga and As nuclei. This creates a subtle 
interplay between electron spin eigenstates affected by the ensemble of nuclear spins (the 
Overhauser shift), nuclear spin states affected by time-averaged electron polarization (the 
Knight shift), and the flip-flop mechanism that trades electron and nuclear spins 

The electron-nuclear interaction has important consequences for quantum information 
processing with confined electron spins yl. Any randomness in the Overhauser shift in- 
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troduces errors in a qubit state, if no correcting measures are taken |5|, |6(. Even worse, 
multiple qubit states, like the entangled states of two coupled electron spins, are redefined by 
different Overhauser fields. Characterization and control of this mechanism will be critical 
both for identifying the problems as well as finding potential solutions. 

Here, the implications of the hyperfine interaction on entangled spin states are studied 
in two coupled quantum dots — an artificial hydrogen molecule — where the molecular 
states can be controlled electrically. A random polarization of nuclear spins creates an 
inhomogeneous effective field that couples molecular singlet and triplet states, and leads 
to new eigenstates that are admixtures of these two. We use transport measurements to 
determine the degree of mixing over a wide range of tunnel coupling, and observe a subtle 
dependence of this mixing on magnetic field. We find that we can controllably suppress the 
mixing by increasing the singlet-triplet splitting. This ability is crucial for reliable two-qubit 
operations such as the SWAP gate 0. 

Furthermore, we find that electron transport itself acts back on the nuclear spins through 
the hyperfine interaction, and time-domain measurements reveal complex, often bistable, 
behavior of the nuclear polarization. Understanding the current-induced nuclear polarization 
is an important step towards electrical control of the nuclear spins. Such control will be 
critical for electrical generation and detection of entangled nuclear spin states 2 > or f° r 
transfer of quantum information between electron and nuclear spin systems 0, 



3] . Even 

more appealing will be reducing the nuclear field fluctuations in order to achieve longer 
electron spin coherence times jiol 11, 3- 

The coupled electron-nuclear system is studied using electrical transport measurements 
through two dots in series jl3j ]. in a regime where the Pauli exclusion principle blocks current 
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flow [14 3- The dots 

are defined using electrostatic gates on a GaAs/AlGaAs heterostruc- 

ture (see Fig. IE and The gate voltages are tuned such that one electron always resides 

in the right dot, and current flows if a second electron tunnels from the left reservoir, through 
the left and right dots, to the right reservoir (see Fig. ID). This cycle can be described using 
the occupations (m, n) of the left and right dots: (0, 1) — > (1, 1) — > (0, 2) — > (0, 1). When an 
electron enters from the left dot, the two-electron system forms either a molecular singlet, 
S(l,l), or a molecular triplet, T(l,l). From S(l,l), the electron in the left dot can move to 
the right dot to form S(0,2). From T(l,l), however, the transition to (0,2) is forbidden by 
spin conservation (T(0,2) is much higher in energy than S(0,2)). Thus, as soon as T(l,l) is 
occupied, further current flow is blocked (Pauli blockade). 

A characteristic measurement of this blockade is shown in Fig. 1A. The suppression of 
current (< 80 fA) in the region defined by dashed lines is a signature of Pauli blockade 



lil Q] (see also Fig. SI and ^|). Fig. IB shows a similar measurement, but with a much 



weaker interdot tunnel coupling t. Strikingly, a large leakage current now appears in the 
Pauli blockaded region, even though the barrier between the two dots is more opaque. 
Furthermore, this leakage current is substantially reduced by an external magnetic field of 
only 100 mT (Fig. 1C). Such a strong field dependence is remarkable because the in-plane 
magnetic field, B ext , couples primarily to spin and the Zeeman energies involved are very 
small (Ez ~ 2.5 /xeV at B ext =100 mT compared to the thermal energy, ~ 15/ieV at 150 
mK, for example). 

Leakage in the Pauli blockade regime occurs when singlet and triplet states are coupled. 
The T(l,l) state that would block current can then transition to the S(l,l) state and the 
blockade is lifted (Fig. ID). As we will show, coupling of singlets and triplets in Figs. 1B,C 
originates from the hyperfine interaction between the electro n sp ins and the Ga and As 
nuclear spins (other leakage mechanisms can be ruled out, see jl7[). 

The hyperfine interaction between an electron with spin S and a nucleus with spin I 
has the form A I ■ S , where A characterizes the coupling strength. An electron coupled to 
an ensemble of n nuclear spins experiences an effective magnetic field ~ Y^^i^-ihi 
with g the electron ^-factor and fiB the Bohr magneton [lj]. For fully polarized nuclear spins 
in GaAs, B^ ~ 5 T For unpolarized nuclear spins, statistical fluctuations give rise 
;o an effective field pointing in a random direction with an average magnitude of 5 T / \fn 
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19]. Quantum dots like those measured here contain n ~ 10 6 nuclei, so \\B^\\ ~ 5 mT. 
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Nuclei in two different dots give rise to effective nuclear fields, B m and B N2 , that are 

> > > 

uncorrelated. Although the difference in field AB N = B m — B N2 is small, corresponding to 

an energy En = <7/x,b|| A.Bjv|| ~ 0.1/xeV, it nevertheless plays a critical role in Pauli blockade. 

The (1,1) triplet state that blocks current flow consists of one electron on each of the two 

dots. When these two electrons are subject to different fields, the triplet is mixed with the 

singlet and Pauli blockade is lifted. For instance, an inhomogeneous field along z causes the 

triplet |T ) = ^(| tl) + I IT)) to evolve into the singlet ^(| Tl> - | IT)). Similarly, the 

other two triplet states, \T + ) = | TT) an d \T~) = | ||), evolve into the singlet due to x and 
> 

y components of ABn- 

The degree of mixing by the inhomogeneous field depends on the singlet-triplet energy 
splitting, Est- Singlet and triplet states that are close together in energy (Est <C En) 
are strongly mixed, while states far apart in energy (E S t E N ) experience only a slight 
perturbation due to the nuclei. 

The singlet-triplet splitting depends on the interdot tunnel coupling t and on the detuning 
of left and right dot potentials Alr- Alr and t are controlled experimentally using gate 
voltages (Fig. IE). V t controls the interdot tunnel coupling. Vl and Vr set the detuning, and 
thereby determine whether transport is inelastic (detuned levels), resonant (aligned levels), 
or blocked by Coulomb blockade (Fig. IF). The coupling of the dots to the leads is held 
constant using V\ ea d- 

The effect of the two tunable parameters t and A LR on the singlet and triplet energies 
is illustrated in Figs. 2A and 2B. For weak tunnel coupling (t ~ 0), and in the absence of a 
hyperfine interaction (En ~ 0), the (1,1) singlet and (1,1) triplet states are nearly degenerate 
(Fig. 2A). A finite interdot tunnel coupling t leads to an anticrossing of S(l, 1) and S(0, 2). 
The level repulsion results in an increased singlet-triplet splitting that is strongly dependent 
on detuning (Fig. 2B). At the resonant condition (A LR = 0, aligned levels), the two new 
singlet eigenstates are equidistant from the triplet state, both with E S t = V%t. For finite 
detuning (but smaller than the single dot S-T splitting), one singlet state comes closer to the 
triplet state (Est ~ t 2 / 'Alr), while the other moves away. Both in Figs. 2A and 2B, singlet 
and triplet states are pure eigenstates (not mixed) and therefore Pauli blockade would be 
complete. 

The additional effect of the inhomogeneous nuclear field is shown in Figs. 2C and 2D. For 
small t (y/2t : t 2 /A LR < E N ), the (1,1) singlet and (1,1) triplet are close together in energy 
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and hence mix strongly (purple lines) over the entire range of detuning. For t such that 
t 2 1 'Alr < En < y/2t, triplet and singlet states mix strongly only for finite detuning. This 
is because Est is larger than En for aligned levels but smaller than En at finite detuning. 
For still larger t (y2t, t 2 /A LR > E N , not shown in Fig. 2), mixing is weak over the entire 
range of detuning. In the cases where mixing between S and T is strong, as in Figs. 2C and 
D (for large detuning), Pauli blockade is lifted and a leakage current results. 

The competition between Est and En can be seen experimentally by comparing ID 
traces of leakage current as a function of detuning over a wide range of t (Fig. 3A). Resonant 
current appears as a peak at Alr = 0; inelastic leakage as the shoulder at Alr > [20]. 
When the interdot tunnel coupling is small, both resonant and inelastic transport are allowed 
due to singlet-triplet mixing, and both rise as the middle barrier becomes more transparent. 
As the tunnel coupling is raised further, the point is reached where Est becomes larger than 
the nuclear field and Pauli blockade suppresses the current (see also Fig. 1A). The maximum 
resonant current occurs at a smaller value of t compared to the maximum inelastic current 
(see inset in Fig. 3A). This is consistent with Est being much smaller for finite detuning 
than for aligned levels, t 2 /A LR < V2t (Figs. 2B, D). 

The experimental knob provided by electrostatic gates is very coarse on the energy scales 
relevant to the hyperfine interaction. However, the external magnetic field can easily be 
controlled with a precision of 0.1 mT, corresponding to a Zeeman splitting (2 neV) that is 
fifty times smaller than En- For this reason, monitoring the field dependence allows a more 
detailed examination of the competing energy scales E S t, E z and E N - 

The competition between Ez and En is clear for small interdot tunnel coupling (Fig. 3B). 
Leakage current is suppressed monotonically with magnetic field, on a scale of ~ 5 mT and 
~ 10 mT for inelastic and resonant transport respectively. The qualitative features of this 
field dependence can be understood from the insets in Fig. 2C. At zero field all states are 
mixed strongly by the inhomogeneous nuclear field, but when E z exceeds En, the mixing 
between the singlet and two of the triplet states (\T + ) and \T_)) is suppressed. An electron 
loaded into either of these blocks further current flow and leakage disappears. 

The magnitude of the fluctuating Overhauser field can be extracted from the inelastic 
peak shape in the limit of small t (as in the inset of Fig. 3B). A fit of the data with a model 



that describes the transport cycle using the density matrix approach 21] is presented in jl7|. 
From this fit, we find a magnitude of the inhomogeneous field y/(AB N 2 ) = 1.73 ± 0.02 mT 
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221. The 



(En = 0.04/ieV), largely independent of Alr over the parameter range studied 
value for the effective nuclear field fluctuations in a single dot is obtained from the relation 
(B N 2 ) = l(AB N 2 ), giving \/ (B N 2 ) = 1.22 mT. This is consistent with the strength of the 
hyperfine interaction in GaAs and the number of nuclei that are expected in each dot jj, . 

The three-way interplay between Est, Ez and En is most clearly visible in the resonant 
current. At an intermediate value of tunnel coupling, t > E N (Fig. 3C), the resonant peak is 
split in magnetic field, with maxima at ±10 mT (see inset). This behavior can be understood 
from the lower inset in Fig. 2D. At B ext = 0, the current is somewhat suppressed compared 
to Fig. 3B, because now Est > En- Increasing B ext enhances the mixing as the \T + ) and 
|T_) states approach the singlet states. The maximum leakage occurs when the states cross, 
at E ST (= y/2t) = E z . Here, E z = 0.25 ± 0.03 ^eV at the current maximum, from which 
we can extract t = 0.18 ± 0.02/ieV for this setting of V t . At still larger B ext , \T + ) and |T_) 
move away from the singlet states again, and the leakage current is suppressed. 

The system enters into a new regime for still higher tunnel coupling (Figs. 3D and 4), 
where it becomes clear that the electron- nuclear system is dynamic. The zero field resonant 
leakage is further suppressed and above 10 mT prominent current spikes appear (left inset). 
The spikes are more dramatically visible in a 3D surface plot of leakage over a broader range 
of field (Fig. 4A). Even for fixed experimental parameters, the current fluctuates in time as 
shown in Fig. 4B. 

We find that time dependent behavior is a consistent feature of resonant transport for 
(Est,Ez) ^> En- For some settings the time dependence is fast (see, e.g., the fluctuations 
of Figs. 4A and 4B), but for others the leakage changes much more slowly. An example 
of the slower time dependence is shown in Fig. 4C. Starting from an equilibrium situtation 
(bias voltage switched off for five minutes), the current is initially very small after the bias is 
turned on. It builds up and then saturates after a time that ranges from less than a second 
to several minutes. This timescale depends on Am, t, and B ext . When no voltage bias is 
applied, the system returns to equilibrium after ~ 80 s at 200 mT. Similar long timescales of 
the nuclear spin-lattice relaxation times have been reported before in GaAs systems and 
quantum dots 2^|. We thus associate these effects with current-induced dynamic nuclear 
polarization and relaxation. 

Evidence that the fast fluctuations too are related to current induced nuclear polarization 
(and cannot be explained by fluctuating background charges alone), is found in their depen- 
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dence on sweep direction and sweep rate 24, 26]. When the magnetic field is swept while 



maintaining fixed A^, the current shows fluctuations at low field but suddenly becomes 
stable at high field (Fig. 4D). The crossover from unstable to stable behavior occurs at a field 
that is hysteretic in sweep direction (Fig. 4D), and this hysteresis becomes more pronounced 
at higher sweep rates (faster than ~ 1 mT/s). The connection between the fluctuations and 
nuclear polarization is also evident from time traces, where instability develops only after 
the nuclear polarization is allowed to build for some time (Fig. S3). 

Unlike the regular oscillations that have been observed in other GaAs structures (see e.g. 



271]). the fluctuations in our measurements are random in time, and in many cases suggest 



bistability with leakage current moving between two stable values. We discuss the origin of 
such fast bistable fluctuations in jl7|. 

The ensemble of random nuclear spins that gives rise to the mixing of two-electron states 



as observed in this experiment also gives rise to an uncertainty of g\iB\/ (Bn 2 ) — 0.03 /ieV in 
the Zeeman energy of one electron. When averaged over a time longer than the correlation 
time of the nuclear spin bath (~ 100 fis) [28], this implies an upper limit on the dephasing 
time of T 2 * = h/ gjiB\l \{Bn 2 ) = 25 ns (following the definition of [^j]), comparable to the 



T 2 * found in recent optical spectroscopy measurements [29]. This value is four orders of 
magnitude shorter than the theoretical prediction for the electron spin To in the absence of 
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nuclei, which is limited only by spin-orbit interactions [30J, 13 1L 13 2j. 

One way to eliminate the uncertainty in Zeeman energy that leads to effective dephasing 



121 ] . Many of the regimes explored in 



is to maintain a well-defined nuclear spin polarization 
this paper show leakage current that is stable when current-induced polarization is allowed 
to settle for some time. These may in fact be examples of specific nuclear polarizations that 
are maintained electrically. 



Supporting online text and figures 

(www. sciencemag.org / cgi / content /full/1113719/DCl) 
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FIG. 1: Pauli blockade and leakage current. (A) Color-scale plot of the current through two 
coupled dots as a function of the left and right dot potentials (voltage bias 800 (leV, Vt = —108 
mV). The experimental signature of Pauli blockade is low current (< 80 fA) in the area denoted 
by dotted grey lines. (B) Analogous data for smaller interdot tunnel coupling (Vt = —181 mV), 
with the same color scale as in (A) . A dramatic increase of leakage current is seen in the lower part 
of the Pauli blockaded area (green/yellow band). Inset: ID trace along the solid grey line, with 
Coulomb blockaded, resonant and inelastic transport regimes marked, see also (F). (C) Analogous 
data for the same tunnel coupling as in (B), but for B ex t = 100 mT. The leakage current from (B) 
is strongly suppressed. (D) Two level diagrams that illustrate Pauli blockade in coupled quantum 
dots (see text). When the (1,1) triplet is changed into the (1,1) singlet (red arrow), Pauli blockade 
is lifted. (E) SEM micrograph showing the device geometry. White arrows indicate current flow 
through the two coupled dots (dotted line). (F) Level diagrams illustrating three transport regimes. 
A: Coulomb blockade; transport would require absorption of energy. □: Resonant transport; the 
dot levels are aligned. +: Inelastic transport; energy must be transferred to the environment, for 
instance by emitting a phonon. These symbols are used to denote inelastic, resonant, and Coulomb 
blockade regimes in Figs. 1, 2, and 3. 
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FIG. 2: Two-electron level diagrams showing energy as a function of detuning Axi?- Detuning 
is denned so that the energy of T(l,l) remains constant as Alr varies (see also Fig. SIB and 
the supporting online text). The panels on the left illustrate the effect of t; the panels on the 
right include the additional effect of an inhomogeneous magnetic field. Pure singlet and triplet 
states are drawn in blue and red respectively, strong admixtures in purple. The blue, white and 
yellow background correspond to the Coulomb blockade, resonant and inelastic transport regime 
respectively. (A) For small tunnel coupling, T(l,l) and S(l,l) are nearly degenerate. (B) For finite 
t, level repulsion between the singlet states results in a larger singlet-triplet splitting compared to 
(A), which depends on detuning. The tunnel coupling does not mix singlet and triplet states. 
For large A^r (but smaller than the single dot S-T splitting) Est ~ t 2 / Alr. (C and D) An 
inhomogeneous field mixes triplet and singlet states that are close in energy (purple lines). For 
clarity only one triplet state is shown in the main panels. (C) For small t, T(l,l) and S(l,l) mix 
strongly over the full range of detuning. (D) For large t, T(l,l) mixes only strongly with the 
singlet for large detuning. The insets to (C) and (D) show the effect of an external magnetic field 
on the two-electron energy levels. All three triplets are shown now; the triplets |T+) and |T_) split 
off from | To) due to B ext . The leakage current is highest in the regions indicated by black dotted 
ellipses (see text). 
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FIG. 3: The measured leakage current results from a competition between En, Est and Ez- 
(A) ID traces of the leakage current as a function of detuning at B ex t = 0, for a wide range of 
tunnel couplings (analogous to the inset of Fig. IB). Coulomb blockade, resonant transport and 
inelastic transport are indicated by colored backgrounds as in Fig. 2. Inset: leakage current along 
the dotted grey and orange lines is shown as a function of Vt- Resonant and inelastic leakage (grey 
and orange markers) reach a maximum at different tunnel couplings {Vt = —190 mV and —150 
mV respectively). (B) For small tunnel coupling (< En), both the resonant and inelastic leakage 
currents drop monotonically with B ex t. Inset: magnetic field dependence of the inelastic current 
along the dotted line (Alr = 40 fieV). (C) For larger t (> En), the resonant leakage current is 
maximum at B ex t =t 10 mT. Inset: field dependence of the resonant peak height (dotted line). (D) 
For still larger t, the resonant current is strongly reduced at low field, then becomes unstable for 
higher field (see inset). 
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FIG. 4: Time dependence of the leakage current reveals the dynamics of the electron- nuclear 
system. This time dependence occurs in the regime corresponding to Fig. 2D. (A) Surface plot of 
electrical transport for Vt = —151 mV. Instability on the resonant peak is visible as sharp current 
spikes. The sweep direction is from + to — Alr, for fields stepped from — to +B ext . (B) Explicit 
time dependence of the resonant current exhibits bistability (Vt = —141 mV, B ex t = 100 mT). (C) 
Lower axis: dynamic nuclear polarization due to electron transport through the device (Vt = —141 
mV, Alr = 0, B ex t = 200 mT), after initialization to zero polarization by waiting for 5 minutes 
with no voltage applied. Top axis: in order to measure the nuclear spin relaxation time, we wait 
for the current to saturate, switch off the bias voltage for a time t re [, and then remeasure the 
leakage current. An exponential fit gives a time constant of 80 ± 40 s (measurements of these long 
timescales result in large error bars, ±20 fA due to 1/f noise). (D) The field dependence of the 
resonant current is hysteretic in sweep direction (Vt = —149 mV). Each trace takes ~ 7 minutes. 



